Stability windows for proto-quark stars 
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We investigate the existence of possible stable strange matter and related stability windows at 
finite temperature for different models that are generally applied to describe quark stars, namely, 
the quark-mass density dependent model, the MIT bag model and the Nambu-Jona-Lasinio model. 
We emphasize that, although the limits for stable strange matter depend on a comparison with the 
ground state of 56 Fe, which is a zero temperature state, the quantity that has to be used in the 
search for strange matter in proto-quark stars is the free energy and we analyze stability windows up 
to temperatures of the order of 40 MeV. The effects of strong magnetic fields on stability windows 
are computed and the resulting mass-radius relations for different stages of the proto-quark star are 
analyzed. 

PACS numbers: 12.39.-x,26.60.-c,95.30.Tg 
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I. INTRODUCTION 

Neutron stars are believed to be the remnants of su- 
pernova explosions They are born hot and rich in 
lcptons. During the very beginning of the evolution, an 
energy of the order of 10 53 ergs is radiated away from 
the neutron star by neutrinos, in a process known as 
deleptonization. Theoretical studies involving different 
possible equations of state that result in different mat- 
ter composition have to be performed. This is because 
the temporal evolution of the star in the so-called Kelvin- 
Helmholtz epoch, during which the remnant compact ob- 
ject changes from a hot and lepton-rich to a cold and 
deleptonized star, depends on two key ingredients: the 
equation of state (EoS) and its associated neutrino opac- 
ity at supranuclear densities 

One class of EOS assumes that the neutron star is 
composed only of hadrons, plus an essential small ad- 
mixture of electrons. In order to construct appropriate 
EOS's for these hadronic stars, one must rely on mod- 
els which describe nuclear matter bulk properties. The 
interior of the neutron star can also include quarks and, 
in this case, two possibilities have been investigated: hy- 
brid stars, constituted by a hadronic phase and a quark 
core, or composed only by charge neutral quark matter 
in ^-equilibrium. In hybrid stars, low density regions are 
composed of hadronic matter, but in high density regions 
deconfinement of the quarks occurs, leading to the forma- 
tion of a quark phase. Many different EOS's have been 
built, including EOS's for the phase immediately after 
the formation of the neutron star, when the neutrinos 
are still trapped, and for the subsequent (deleptonized) 
phases after the neutrinos escape. 

The assumption underlying the existence of quark stars 
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is based on the Bodmer-Witten conjecture @, Q ■ These 
authors have claimed that it is possible that the inte- 
rior of a neutron-like star does not consist primarily 
of hadrons, but rather of strange matter (SM), which 
is composed of deconfined quarks, including up, down 
and strange quarks, plus the leptons necessary to en- 
sure charge neutrality and /3-equilibrium. This possibility 
arises because a phase transition from hadronic to quark 
phase is possible at densities present in the interior of 
neutron stars. It has been argued @ that strange matter 
is the true ground state of all matter. If this is the case, as 
soon as the core of the star converts to the quark phase, 
the entire star converts. SM was first considered in cal- 
culations obtained within the MIT bag model framework 
[13, [nj. More sophisticated treatments for SM, based on 
the Nanbu-Jona-Lasinio fl2l - [l6j . the color flavor locked 
phase [l7l - [2l| and the quark-mass density dependent [22[ 
models also exist in the literature. 

An important ingredient in the SM hypothesis is the 
stability window, identified with the model parameters 
that are consistent with the fact that two-flavor quark 
matter (2QM) must be unstable (i.e., at zero tempera- 
ture its energy per baryon has to be larger than 930 MeV, 
the iron binding energy) and SM (three-flavor quark mat- 
ter) must be stable (i.e., its energy per baryon must be 
lower than 930 MeV, also at T = 0). This problem has 
already been discussed in the literature, but we believe 
that some aspects deserve more investigation and this is 
the main goal of the present work. First of all, at zero 
temperature, the stability window is normally obtained 
for neutral matter in /3-equilibrium, as for instance in 
Ref. [23l - [25l ] . Nevertheless, one has to bear in mind that 
stable nuclear matter (as in iron) is not charge neutral 
and does not contain electrons. Actually, its proton frac- 
tion is Y p = 0.46, very close to symmetric matter and 
this is a good reason to analyze also matter with equal 
quark chemical potentials, as done in Ref. [2(| for zero 
temperature systems. Moreover, an important charac- 
teristic of the deconfinement transition is that deconfined 
quark matter is transitorily out of equilibrium with re- 
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spect to weak interactions [20, 27, 28]. It was also shown 
in Ref. [26| that slightly different results are obtained if 
one considers matter with identical quark chemical po- 
tentials, corresponding to equal quantities of u and d 
quarks in two flavor-matter, or charge neutral matter in 
/3-equilibrium, as expected in stars. For the models un- 
der investigation, the stable parameter region is larger 
if the /3-equilibrium condition and charge neutrality are 
imposed. 

Another point worth mentioning refers to finite size ef- 
fects, very often disregarded in the literature. As shown 
in Ref. if finite size effects are taken into account, 
iron binding energy is around 4 McV lower, i.e. 934 McV. 
However, one can see (2(| that even if finite size effects are 
taken into account, the stability window remains practi- 
cally unchanged. 

Understanding the stability of strange matter in proto- 
quark stars is also important because it is known that 
stars cool down slowly after deleptonization. According 
to numerical simulations, during the first tens of seconds 
of evolution the proto-ncutron star cools from T = 40 
McV (in the center) to temperatures below 2-4 MeV . 
One does not expect that proto-quark stars are constant 
in temperature. Actually, it is believed that the entropy 
per baryon is fixed, so that the temperature varies in the 
interior of the star, according to its density. A reasonable 
assumption is that during the evolution process, differ- 
ent snapshots can be simulated though different entropies 
per baryon and trapped neutrinos. See Refs. [29|, H(| for 
some examples of studies of isolated star evolution. This 
analysis is important as it has been shown that evolution 
of isolated stars differs from the one of stars in binary 
systems throughout the whole star life [HI, |32| . 

The stability window at finite temperature has already 
been discussed in the past f33T - l36j . In all these papers, 
except in Ref. [33[, the importance of the free energy 
was disregarded, but even in Ref. [33^ . where it was in- 
troduced in the calculations, the quantity used in the 
search for stable strange matter was always the binding 
energy per baryon. Of course, as expected from calcula- 
tions in the macrocanonical or grand-canonical ensemble, 
the quantity related to the thcrmodynamical potential is 
the free energy per baryon (T/A = f / p B = (e— Ts)/ p B ), 
where / is the free energy density, pb the baryon density, 
e the energy density, T the temperature and s the entropy 
density of the system [37| • For zero temperature sys- 
tems, the free energy density becomes the energy density 
and hence, the binding energy per baryon {B/A = e/p B ) 
is analyzed at zero temperature in the search for stable 
matter. The choice of appropriate parameters compatible 
with stable SM at finite temperature systems requires, in 
addition to the investigation performed at zero tempera- 
ture, a careful study of the free energy per baryon. A very 
similar situation is seen when one looks for the pos sible 
existence of the pasta phase in finite systems [38ll39| . Of 
course, the lower limits of the stability windows, given 
by the points where the two-flavor quark matter has an 
energy per baryon larger than 930 McV, remain the same 



because the 56 Fe ground state is a zero temperature sys- 
tem. 

Another important aspect related to neutron stars is 
their magnetic fields. While magnetars are expected to 
bear extremely high m agnetic fields of up to B = 10 15 G 
on the surface 



some strength 



-|43j , all pulsars have magnetic fields of 
44| . For this reason, any complete analy- 
sis of pulsar features, including stability windows, should 
take into account magnetic field effects. Note that when- 
ever tackling this problem, anisotro py is a matter that 
has to be carefully considered [2l|, l30l. l45l - [55j , once it may 
be important for ultra high magnetic fields. In this work, 
however, we are not going to take into account anisotropy 
of matter, since we are not primarily interested in ana- 
lyzing macroscopic properties of stars, but general trends 
of different microscopic models in specific situations. 

For the above mentioned reasons, it is appropriate to 
have a better understanding of stability windows also at 
finite temperature (including the case of fixed entropy 
per baryon) subject to magnetic fields. This is the scope 
of the present paper, organized as follows: in Section 
II, we briefly discuss the formalism and the EOS's of 
different models; in Section III we present the results 
and conclusions and in the last section, the final remarks 
are drawn. 



II. FORMALISM 

We discuss next the following models for proto- 
quark stars: the quark-mass density dependent (QMDD) 
model, the MIT bag model and the Nanbu-Jona-Lasinio 
(NJL) model. We restrict ourselves to the main formulae 
and refer the reader to previous references for detailed 
calculations. 



A. QMDD model 

We start from the QMDD model [H, [H, [H HHH, 
which is based on a phcnomcnological approach where 
the dynamical masses of the three lightest quarks scale 
inversely with the baryon number density 
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where C is the constant energy density in the zero quark 
density limit. 

The pressure, energy density and baryonic density of 
the system are respectively given by 
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and 7, is the degeneracy of each quark i = u,d, s 
taking into account spin and number of colors. 
The distribution function for quarks (/+) and anti- 
quarks (/_) are the Fermi-Dirac distributions f±i = 
[1 + exp [(E* =p Hi) /T]] -1 where ^ is the chemical po- 
tential of each particle species and E* (p) = yjk 2 + to* 2 . 

This model presents a thermodynamical inconsistency 
which was tackled in different ways in the literature and 
resulted in different versions of the QMDD model. This 
problem was discussed in many papers J2(| l57l - [59| and we 
restrict ourselves to two versions: 1) the most commonly 
used [23| , where the pressure at the density corresponding 
to the minimum of the free energy per baryon could be 
non-zero, depending on the matter studied (SM or 2QM) 
and to which we refer next as version 1 (QMDDvl) and 
2) the one that presents a remedy to the thermodynami- 
cal inconsistency, in such a way that the minimum of the 
ener gy p er baryon corresponds to the point of zero pres- 
sure [58|, and to which we refer as version 2 (QMDDv2). 
The confinement parameter values for which the stability 
windows are obtained in different versions of the model 
at zero temperature are different, as shown in Ref. f26j . 



B. MIT bag model 

For the MIT bag model [l(| E| the quark masses are 
fixed. Therefore, expressions @, (JS|) and ((H) have m* 
replaced by m and B{ns, f^i) by a bag constant B. We 
take advantage of the simplicity of this model and include 
magnetic field effects in the calculation. In this case, 
the main expressions are written next (see Ref. [3(| for 
details). For the pressure, energy density and baryonic 
density they read, respectively 



=(/+* + /-i)-B, 



(6) 



EE ^2 \H \ B J dky/kf + mSiUi + f^ + B, 



^=Ef = JEE^W dhU +i -f-i), (8) 

i i v 

where 7,; is once more the degeneracy of each quark 
i = u,d,s taking into account spin and number of col- 
ors. The difference now is that both spin projections 



contribute for Landau levels v > 0, but only one of them 
contributes for v = 0. Note that we only consider the case 
where the magnetic field is parallel to the z direction so 
the energy levels in the x and y directions are quantized. 
In this case qi is the electric charg e of each quark , B is the 
magnetic field strength, rftj = y/rn? + 2\qi\Bv the mass 
of each quark {m Ut d — 5 McV, m s can take different val- 
ues) modified by the magnetic field and Ei = \/k 2 + rhi 2 . 



C. NJL model 

For the NJL model (see Ref. [l2l . [l3j for details) , we 
need to evaluate the grand-canonical thermodynamical 
potential for the three flavor quark sector 



Q= -P = e-Ts- 



(9) 



where ensures that Q = in the vacuum. As a re- 
sult, we obtain for the pressure and entropy density of 
the system in the mean field approximation in the pres- 
ence of a magnetic field the expressions given next. See 
Appendix for the formulas without magnetic field effects 
and refer to Refs. (60l - [62j for details on the calculations 
of the expressions given below: 
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(11) 

The energy density of the system can be calculated from 
Eq. (J5J taking the baryon density to be the same as in 
Eq.([8j). For this model we split the degeneracy of each 
quark into the spin degeneracy a v (again both projec- 
tions contribute for v > but only one for v = 0) and 
color degeneracy N c . In the above equations, G and K 
are coupling constants, qi is the electric charge of each 
quark, B is the magnetic field strength and the con- 
tribution from the gas of quasiparticles for each flavor 
0i = (er c + 0" lag + e T ed ) Mi contains 3 different contri- 
butions: the vacuum, the magnetic and the medium one 
given by 
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FIG. 1: (Color online) Stability windows of SM obtained with 
the MIT model through the analysis of the free energy per 
baryon and the binding energy per baryon (both ~ 930 MeV) 
shown for different temperatures. 



FIG. 2: (Color online) Stability windows of SM obtained with 
the QMDDvl model through the analysis of the free energy 
per baryon and the binding energy per baryon (both ~ 930 
MeV) shown for different temperatures. 



Above, we have defined ea = \/A 2 + Mf with A rep- 
resenting a non covariant ultra violet cut off |63j . Xi = 
Mf/{2\ qi \B) and C (1 ' 0) (-l,x 4 ) = d({z, Xi )/dz\ z= ^ with 
(^(z,Xi) being the Ricmann-Hurwitz zcta function. 

Each of the quark condensates, 4>i = (q~iqi) = (</>™ c + 
^mag _|_ ^medj^ a i so corr k ams 3 different contributions: 

the vacuum, the magnetic and the medium one given by 
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where E* 



y/k'f + s t (v,B) 2 and Si{v,B) 



^2\qi\Bv is the constituent mass of each quark 
modified by the magnetic field. Once more we consider 
the case where the magnetic field is parallel to the z di- 
rection so the energy levels in the x and y directions 
are quantized as done in Refs. (60l - [62^ . Minimizing the 
grand-canonical thermodynamical potential fi with re- 
spect to Mi leads to three gap equations 

Mi = m, -4G0.+ 2K^ k , (18) 

with cyclic permutations of i, j, k. 

D. Chemical equilibrium 

In the description of compact stars, both charge neu- 
trality and chemical equilibrium conditions have to be 



imposed [H, |64| . The first condition can be written for 
quarks and lcptons as 

2p« = Pd + Ps + 3 (p e + Pti) , (19) 
and the second condition can be written as 

(l s = /Lt d = M« + Me, Me = /V, (20) 

where the lepton densities can be calculated through 
Eq. ([4]) or Eq. ([8]), depending on the inclusion or not 
of the magnetic field, with appropriate substitutions for 
the masses and electric charge. The lepton masses are 
m e = 0.511 MeV and m M = 105.66 MeV. For the elec- 
tron and muon pressure and energy density we use Eqs. 
(0), and © or © and ([7]), in the case without and with 
the inclusion of the magnetic field, respectively, all with 
B = 0. 

In earlier stages, when the neutrinos are still trapped 
in the interior of the star, Eq. (|20p is replaced by 



Ms = Md = M« + Me 



1-1,, 



Mm = Me 



and 



= M^ • 



(21) 



(22) 



The independent chemical potential [i ve appears as the 
system gains one conserved quantity, the lepton fraction 
Yl = (Pe + Pu + Pv„ +Pv„)/pb, which, according to simu- 
lations [H, [65j], can reach Y\ = 0.4. Note that the degen- 
eracy factor for the neutrinos is 7; = 1. 



III. RESULTS AND CONCLUSIONS 

We start by analyzing the stability window related to 
proto-quark stars described by SM obtained with the 
MIT bag model. As stated in the Introduction, instead 
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FIG. 3: (Color online) Stability windows of SM obtained with 
the QMDDv2 model through the analysis of the free energy 
per baryon shown for different temperatures. 



Set 


Model 


m a (MeV) 


S 1/4 (MeV) 


C (MeV/m -3 ) 


A 


MIT 


150 


155 


N/A 


B 


QMDDvl 


150 


N/A 


78 


C 


QMDDvl 


100 


N/A 


85 


D 


QMDDv2 


150 


N/A 


194 


E 


QMDDv2 


100 


N/A 


215 



TABLE I: Parameters for the MIT and QMDD models. 



of considering the binding energy, the quantity that has 
to be studied in obtaining the upper limit of the stability 
window is the free energy. In order to show how differ- 
ent the results are, we plot both the binding energy per 
baryon and the free energy per baryon in Fig Q] as a func- 
tion of the bag constant and the strange quark mass for 
the MIT bag model, without the inclusion of magnetic 
field effects. Notice that we have used, for two-flavor 
quark matter (2QM), the fact that fj, u = Pd, which gives 
symmetric matter (p u = pd) and, to be consistent, for 
SM we have used p u = [id = /j, s . As the strange quark 
mass is much larger than the masses of quarks u and d, 
its relative density is considerably lower. For T = 0, the 
binding energy per baryon has been calculated numer- 
ically for 2QM and SM, respectively, and SM is stable 
in the shaded region shown in Fig[TJ The lower limit, 
vertical straight line, is due to the requirement that two 
flavor quark matter is not absolutely stable. We have 
repeated the calculation for temperatures up to 40 MeV, 
as they are relevant for the proto-neutron star evolution 
simulations. Our results reproduce the same behavior as 
the ones given in Refs. [H, [35|,[36| for the binding energy 
per baryon, i.e., as temperature increases, the stability 
windows move to the left, i.e., becomes more restrictive. 
However, when we consider the free energy per baryon, 
the stability window moves to the right and a wider range 
of constants becomes possible to ensure stable matter, 
contrary to what has been published so far. 
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FIG. 4: (Color online) Stability windows of SM obtained with 
the MIT model with parametrization A without (top) and 
with (bottom) the inclusion of a magnetic field of 7.2 x 10 18 
G shown for different temperatures. The orange rectangle 
shows the highest temperature that overlaps with T = 0. 

In order to choose adequate values for the strange 
quark mass and the constant that enters in the quark 
masses in the QMDDvl model, we display stability win- 
dows for finite temperatures in Fig. [2] The shaded 
area corresponds once more to the zero temperature case. 
Note that considering the free energy per baryon as the 
criterion for stability, the windows move right with the in- 
crease of the temperature, as in the MIT model. We anal- 
yse also the QMDDv2 model and the results are shown 
in Fig. [3] One can see that the values of the confinement 
constants are somewhat different for this version of the 
model, with obvious consequences on the resulting stellar 
matter configuration, as will be seen next. 

We then fix the strange quark mass to 150 MeV 
(parametrization A in Table HJ and calculate the over- 
lap of stability windows at various temperatures through 
the analysis of the free energy per baryon for the bag 
model. In Fig. @] we display results without and with the 
inclusion of the magnetic field. For each temperature, 
we have used the prescription given above for unstable 
2QM and stable SM to restrict the bag model parameter 
values. We can see that after a certain temperature, i.e., 
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FIG. 5: (Color online) Stability windows of SM obtained with 
the QMDDvl model with parametrizations B (top) and C 
(bottom) without the inclusion of magnetic field effects shown 
for different temperatures. The orange rectangle shows the 
highest temperature that overlaps with the zero temperature 
one. 



T = 33.4 MeV for non-magnetized matter and T = 21.7 
MeV for magnetized matter, the boxes representing the 
parameter values stop overlapping. However, consider- 
ing that the lowest limit should always be established for 
the T = 0/B = case, a wide range of parameters can 
be used for finite temperature matter inside stars. For 
example, if we consider matter with fixed entropy per 
baryon S/A = 2, the temperature can reach 40 MeV in 
the center of the proto-quark star. In this case, £> 1//4 can 
vary from 147 to 168.7 MeV, for unmagnetized stars and 
147 < B 1/4 < 172.2 MeV for magnetized stars. 

In Fig. [5] the same calculation is performed for the 
QMDDvl model using parametrizations B and C from 
Table HI In the first case, one can see that for T = 30 
MeV, no overlapping would be possible within this model 
if the lower limit of the stability window were obtained 
also considering the finite temperature system. This hap- 
pens because there is an overall smaller overlap of stabil- 
ity windows for finite temperature with the QMDD than 
with the MIT bag model. The last temperature where 
some overlap still exists with the zero temperature box 
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FIG. 6: (Color online) Stability windows of SM obtained with 
the QMDDv2 model with parametrizations D (top) and E 
(bottom) without the inclusion of magnetic field effects shown 
for different temperatures. The orange rectangle shows the 
highest temperature that overlaps with the zero temperature 
one. 



is obtained for T = 17.8 MeV for parametrization B and 
T = 31.6 MeV for parametrization C. 

For the sake of completeness, we recalculate the stabil- 
ity windows for the QMDDv2 model and the results are 
shown in Fig. [5] The interpretation of the appropriate 
parameter values for stellar matter follow the explana- 
tion given above. Notice that, as expected from the re- 
sults displayed in Fig. [5J the overall parameter values are 
very different, making the equations of state of QMDDv2 
much softer than the ones of QMDDvl. Comparing Figs. 
[5] and [51 one can see that the maximum values allowed for 
the parameter C within the QMDDvl for T = 40 MeV 
are much lower than the minimum values allowed for the 
QMDDv2 model at zero temperature. 

From hereon all calculations for the MIT bag model 
and the QMDD model are made considering the sets of 
parameters shown in Table [J and discussed above. These 
parametrizations have realistic strange quark masses and 
they are the ones that comprehend stability windows for 
the largest possible range of temperatures. 
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Set 


Model 


S/A 




r CIO 18 G) 


M„,„„ (MrO 
lyJ -max V lvA y ) 


R (km) 




T„ CMeV) 


A 


MIT 


n 


N/F 


n 

u 


1.62 


Q 01 


8 9^ 


u 


A 

I i 


MIT 


I 


i 


n 

V.J 


1.64 


Q 1 


7 Qfi 


-LO.OZ 


A 


MIT 


2 


0.4 


o 


1.65 


9.15 


7.85 


26.59 


A 


MIT 


o 


N/F 


6.64 


2.02 


9.04 


8.31 


o 


A 


MIT 


1 


0.4 


4.71 


1.95 


9.05 


8.82 


12.20 


A 


MIT 


2 


0.4 


4.44 


1.93 


9.08 


8.69 


24.27 


B 




n 

u 


N/F 


n 

u 


2.28 


1 9 0^ 




n 


B 


oivinnvi 


I 


0.4 


o 


2.31 


12.16 


4.39 


11.34 


B 


QMDDvl 


2 


0.4 


o 


2.33 


12.19 


4.46 


22.80 


c 


QMDDvl 





N/F 


o 


2.26 


11.76 


4.65 





c 


QMDDvl 


1 


0.4 


o 


2.28 


11.75 


4.76 


11.78 


c 


QMDDvl 


2 


0.4 


o 


2.29 


11.76 


4.79 


23.62 




OMDDv2 


o 


N/F 


o 


1.60 


8.42 


9.37 


o 


13 




I 


4 


n 

u 


1.62 


8 dp, 




1 1 SO 


D 


QMDDv2 


2 


0.4 





1.62 


8.46 


9.23 


27.64 


E 


QMDDv2 





N/F 





1.59 


8.22 


9.75 





E 


QMDDv2 


1 


0.4 





1.58 


8.16 


9.93 


14.24 


E 


QMDDv2 


2 


0.4 





1.58 


8.16 


9.92 


28.51 



TABLE II: Output results given by the numerical solution of Tolman Oppenheimer Volkoff (TOV) equations for the MIT and 
QMDD models at different snapshots of the proto-quark star evolution. For the zero temperature case, the lepton fraction is 
not fixed (N/F). 



Finally, we have used EOS's calculated with the pa- 
rameters shown in Table U as input to the TOV equations 
[6ol |67| . The resulting mass radius relations are plotted 
in Fig. [7] The first model analyzed is the MIT at dif- 
ferent entropies per baryon, indicating different evolution 
stages of the star without and with magnetic field effects. 
In the cases in which the magnetic field is taken into ac- 
count, it is considered to be variable, i.e., its magnitude 
increases with baryon chemical potential (density) fol- 
lowing Refs. j68l - l7l| ) from a surface value lower of 10 15 
G to a value lower than 7.2 x 10 18 G in the center of the 
star. Additionally, for this specific calculation, we add 
the pure magnetic field contribution proportional to B 2 
to the pressure and energy density. As seen in many pa- 
pers (as an example, see Refs. [3(1 EH), the magnetic field 
makes the EOS stiffer (when the pressure perpendicular 
to the magnetic field is used as input) and consequently, 
more massive stars are obtained in this case. This re- 
sult is consistent with results from realistic simulations 
including magnetic fields [7^, [74[. Note that when we 
include magnetic fields, the effect of a higher entropy per 
baryon (and consequently temperature) is to lower the 
maximum mass a star can have, instead of increasing it 
due to thermal effects. As pointed out in Ref. j3(j, this is 
due to change in baryon number between different evo- 
lution stages, only allowed in stars belonging to binary 
systems. 

For the QMDDvl and QMDDv2 models, the results 
for different parametrizations are also shown in Fig. [7] 
While QMDDvl always reproduces more massive stars 



than the MIT bag model due to the extra terms appear- 
ing in the grand-canonical thermodynamical potential, 
QMDDv2 models present results that are similar to the 
MIT ones, since the EOSs are much softer than the ones 
for QMDDvl. One can see that the QMDDvl model can 
reproduce very massive stars, as the ones recently mea- 
sured [iBI, [7(| , but fails to describe stars with low radii, 
if their existence is confirmed. On the other hand, the 
QMDDv2 model can describe stars with radii even lower 
than the MIT model. For different values of the con- 
finement constant, it was already shown in Ref. [2(| that 
this version of the QMDDvl model generates larger max- 
imum masses thant the MIT model systematically, for all 
possible parameters inside the stability window. These 
results are not changed even when strong magnetic fields 
are considered in the MIT model. 

Furthermore, we include fixed lepton fraction in the 
EOS's in order to improve our description of proto-quark 
stars. For the QMDD models, only electrons and their 
related neutrinos are considered within the lepton frac- 
tion, while in the MIT model, muons and their neutrinos 
are also taken into account. The results for the QMDD 
and MIT bag model are shown in Table [TT1 For the later, 
we also show the inclusion of (variable) magnetic field ef- 
fects. Note that the central magnetic field reached within 
the MIT bag model is higher at later stages of the star 
evolution, when the stars are also more massive due to a 
higher baryon number. In this table, besides maximum 
mass and respective radius, we also show central energy 
density and central temperature for each case considered. 
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FIG. 7: (Color online) Mass-radius relation for different mod- 
els shown for different entropies per baryon. We have also 
included lower and upper limits of the masses and radii of 
EXO 0748-676 and 4U 1608-52. The shaded clouds refer to 
the la and 2a confidence ellipse of the results obtained in for 
the EXO 1745-248 \?Mnk- The four different sets of curves 
are presented in the legend going from higher star masses to 
lower star masses. 



We now analyze the NJL model, for which the param- 
eters are fixed so as to fit the values in vacuum for the 
pion mass, the pion decay constant, the kaon mass and 
the quark condensates. We consider the sets of parame- 
ters given in Tabic Hill taken from Ref. [8(| and references 
therein. In Fig. HI the results for the SU(2) version of the 
model are shown. All four parametrizations are consid- 
ered for T = and non-magnetized matter. We notice 
that parametrization 1 results in an unbound system. 
The other three parametrizations yield bound systems, 
but the binding energies are larger than the one for iron. 
Therefore, the SU(2) version of the NJL model cannot 
be stable without the inclusion of the strange quark, as 
expected. Next, we include the strange quark in the cal- 
culation. 

The results for both SU(3) parametrizations 5 and 6, 
shown in Table Hill are seen in Figs. 151 andl 101 respectively. 
A detailed discussion on the possibility of stable strange 
quark matter with the NJL model at zero temperature 
was done in section 3.3.2 of Ref. [80j . where the author 
claims that the NJL model does not support the idea of 
absolutely stable quark matter. However, in the present 
work we can see that as the temperature increases, the 
free energy per baryon decreases, but the system becomes 
unbound. On the other hand, if we consider the influence 
of the magnetic field, the system becomes more bound 
as the magnetic field strength increases and, therefore, 
more stable. The latter had already been pointed out in 
Refs. (8ll - l83l | for zero temperature. The dots show the 
points where the pressure is zero. These points do not 
coincide with the minima of the free energy per baryon 




FIG. 8: (Color online) Binding energy per baryon as a func- 
tion of baryon density shown for different parametrizations 
of the SU(2) NJL model at zero temperature. The dots are 
points of zero pressure. 



Set 


Group 


A [MeV] 


GA 2 


m u , d [MeV] 


KA 5 


m s [MeV] 


1 


SU(2) 


664.3 


2.06 


5 


N/A 


N/A 


2 


SU(2) 


587.9 


2.44 


5.6 


N/A 


N/A 


3 


SU(2) 


569.3 


2.81 


5.5 


N/A 


N/A 


4 


SU(2) 


568.6 


3.17 


5.1 


N/A 


N/A 


5 


SU(3) 


602.3 


1.835 


5.5 


12.36 


140.7 


6 


SU(3) 


631.4 


1.835 


5.5 


9.29 


135.7 



TABLE III: Parameters for the NJL model obtained from Ref. 
18011 and references therein. 



for large magnetic fields. 

Note that for strong magnetic fields and low tempera- 
tures we see metastable configurations in Figs. 151 and 1TU1 
Those are related to the Haas-van Alphcn oscillations 
[8~4T - l87j that usually appear within these limits. These 
oscillations are related to when the Fermi energies of the 
charged particles cross the discrete threshold of a Landau 
level. See Ref. [§H for more details on the effect of mag- 
netic field on the Haas- van Alphen oscillations within the 
NJL model. 



IV. FINAL REMARKS 

We presented in this paper a review of stability win- 
dows for quark matter at finite temperature using differ- 
ent models. Wc pointed out that the correct quantity to 
analyze in order to obtain the upper limits of the stabil- 
ity windows is the free energy per baryon, instead of the 
binding energy per baryon. In this case, for the QMDD 
as well as for the bag model, the correct values of the re- 
spective coupling constants (C and S 1 / 4 ) assume larger 
values at larger temperatures. Even when a strong mag- 
netic field is applied to the system, this behavior does 
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FIG. 9: (Color online) Free energy per baryon as a function 
of baryon density for parametrization set 5 of the SU(3) NJL 
model shown for different temperatures. Full lines represent 
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pressure. 



1100 
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2 x 10 G. The dots are points of zero 
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FIG. 10: (Color online) Free energy per baryon as a function 
of baryon density for parametrization set 6 of the SU(3) NJL 
model shown for different temperatures. Full lines represent 
calculations with B — 1 x 10 18 G, dashed lines calculations 
with B = 5 x 10 18 G, dotted lines B = 1 x 10 19 G and dot- 
dashed lines B = 2 x 10 19 G. The dots are points of zero 
pressure. 



not change, but it is enhanced. 

Based on the results stated above, we chose vari- 
ous parametrizations for the the QMDD model and two 
paramctrizations for the bag model. For the latter they 
correspond to the cases without and with magnetic field 
effects. With these EOS's in hand, we simulate proto- 
quark stars for fixed entropies per baryon and fixed lep- 
ton fraction. We see that in all cases analyzed (with- 
out and with magnetic field effects), the QMDDvl re- 
produces more massive stars than the MIT bag model, 
while QMDDv2 produces stars with low masses and low 



radii. 

The NJL model allows a more realistic description of 
quark matter, as it contains chiral symmetry restora- 
tion/breaking. For these reason, we studied both versions 
of the model, SU(2) and SU(3) in detail to check whether 
the stability for the existence of SM is satisfied with the 
parametrizations normally used in the literature. The 
first group is not stable, as it is not bound bellow 930 
McV. The second group has stable paramctrizations for 
different temperatures (found by the analysis of the free 
energy per baryon), specially when strong magnetic fields 
are considered. 
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Appendix 

Here we present the most relevant formulas for the 
NJL model without the inclusion of magnetic field ef- 
fects. Once more, we start from Eq. (j9)), from which we 
obtain for the energy density and entropy density of the 
system in the mean field approximation 



-2N c y / - — 

c 4-y (2tt 



d 3 ki kf + m i M l 



E, 



(f-i - /+*MA 2 - kf) 



2G(<pl + 4> 2 d + cj>f) + IK^Ms - eo, 



(23) 



= -2^E 



d 3 k, 



(27T) 



f+i In (/+,,) + (1 - /+,) 



x In (1 - f +i ) + f_t In (/_<) + (1 - f_i) In (1 - /_.,) 
x 6(A 2 -kf), (24) 



where E* = -J kf + Mf and Mj is the constituent mass 
of each quark. Minimizing the grand-canonical thcrmo- 
dynamical potential Vl with respect to Mi leads to the 
three gap equations Eq. (fT5)l . The baryon density has 
the same form as given in Eq. @. 

The quark condensates are defined, for each of the fla- 
vors i = u, d, s as 



4>% = (n <n) 



-2N c y / - — 

(2tt 



d 3 h Mi 



(2tt)3 Ei 
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